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Abstract. Spin-1/2 magnets with kagome geometry, being for years a generic object 
yJ , of theoretical investigations, have few real material realizations. Recently, a DFT- 

;_( ■ based microscopic model for two such materials, kapellasite Cu3Zn(OH)6Cl2 and haydeeite 

TJ^ ' Cu3Mg(OH)6Cl2, was presented [Janson O, Richter J and Rosner H, arXiv:0806.1592 . Here, 

we focus on the intrinsic properties of real spin-1/2 kagome materials having influence on the 

C^ ' magnetic ground state and the low-temperature excitations. We find that the values of exchange 

C I integrals are strongly dependent on O — H distance inside the hydroxyl groups, present in most 

I ■ spin-1/2 kagome compounds up to date. Besides the original kagome model, considering only the 

' O ' nearest neighbour exchange, we emphasize the crucial role of the exchange along the diagonals 

C of the kagome lattice. 

o 

Two-dimensional (2D) magnets with a kagome lattice arrangement of magnetic ions (Figured]) 
K^ ' are geometrically frustrated due to a triangular-like arrangement of nearest neighbours (NN) 

Jifl^ ■ leading to an unusual highly degenerate classical ground state (GS) [T]. Strong quantum 

(<— ^ I fluctuations arising from spin-1/2 ions can lift the classical degeneracy and drive the system 

psj ■ into a magnetically disordered quantum paramagnetic state, which might be applied for future 

[^ . quantum computational applications [2]. Nonetheless, in real materials couplings to further 

^^ I neighbours [31 \^ and between the kagome layers [5] are always present and influence the GS and 

the thermodynamics. 

In a recent paper [^ we have performed DFT calculations for two new natural isostructural 
materials with spin-1/2 kagome lattice — kapellasite Cu3Zn(OH)6Cl2 |6j and haydeeite 
S^ . Cu3Mg(OH)6Cl2 O [8]. The local density approximation (LDA) yields a metallic solution 

J^ I (left panel in Figure [2|) in contrast to the experimentally observed insulating behaviour due 

to underestimation of strong on-site correlations for the Cu^"*" sites. To account for this deficit, 
we carried out an effective one-band tight-binding (TB) fit of the LDA bands (Figure [21 right 
panel) and subsequently mapped the resulting transfer integrals tj to an extended Hubbard 
model and to a Heisenberg model with Jj = Atj/Ueff- This mapping provides an estimate for 
the antiferromagnetic (AFM) part of the exchange. To obtain the values for the total exchange, 
consisting of AFM and ferromagnetic contributions, we have performed supercell calculations 
for various spin arrangements. The correlations were treated in the LSDA-|-f/ approximation. 
We have found that besides the standard kagome model based on NN interactions (Ji), an 
additional coupling along diagonals of kagome hexagons (J^) is relevant (Figure [1]). Since the 
absolute values of Ji and Jd change only the overall temperature scale, their ratio a=Jd/ Ji is 
crucial for the GS and low-energy excitations. 
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Figure 1. Kagome layers in 
the crystal structure of kapel- 
lasite (haydeeite) . CUO4 pla- 
quettes are shown in yellow, 
ZnOg (MgOg) octahedra are 
shown in blue. The rele- 
vant exchange paths are high- 
lighted. 
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Figure 2. Left panel: total and atom-resolved densities of states of kapellasite Cu3Zn(OH)6Cl2 
(left) and the corresponding band structure (right). Right panel: band structure of antibonding 
dpa bands crossing the Fermi level and the effective one-band tight-binding fit for different O — H 
distances (experimental: 0.78 A, LDA bands — dashed lines, TB fit — green squares; optimized: 
1 A, LDA bands — solid lines, TB fit — red circles). 



One focus of this paper is to emphasize the influence of side groups (which are commonly 
neglected in model physics) on different exchange paths. A key issue of theoretical investigations 
of real materials is the construction of a relevant model, which describes experimental data and 
has predictive power. For cuprates, that are insulators with a 3d^ Cu^+ configuration, the 
extension from the initial Heisenberg idea to consider only spin degrees of freedom (Heisenberg 
model) to a multiband Hubbard model with a separate treatment of Cu ?>d and O 2p orbitals 
became computationally feasible only recently. On the other hand, the famous Goodenough- 
Kanamori- Anderson (GKA) rules formulated quite early on an empirical basis provide a simple 
intuitive picture for the exchange interactions based on geometrical quantities. Due to the 
simplicity of the GKA rules, they are violated in many cases. A natural way to search for a 
possible origin of these phenomena is to improve the model by including ligand fields effect. In 
general, such models are very difficult to evaluate. Even in case of simple systems like CuGeOa 
the models are quite complicated and contain many parameters [9]. 

Here, we present the results of a DFT-based modeling for Cu3Zn[Mg](OH)6Cl2. In these 
compounds, the leading exchange interactions are strongly dependent on the position of H 



Figure 3. Finite kagome lattice of A^ = 36 sites. The solid lines represent the NN bonds Ji 
and the dashed lines the diagonal bonds J^. The coloured circles indicate the spin orientations 
of the twelve-sublattice classical GS relevant for J^i > 0. Along the chains formed by diagonal 
bonds there is an antiparallel (Neel) spin alignment (e.g. the spins at the light and dark green 
circles). On each triangle formed by NN bonds Ji (e.g. spins on sites 21, 22, 23) there is a 120° 
spin arrangement. In addition, two antiparallel (Neel) spin-sublattices are perpendicular to one 
other group of two Neel-like sublattices, i.e. •J.^ ;^_L^ ;^_L^ 
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Figure 4. GS spin-spin correlation 
(SqSj) in dependence on a = Jd/Ji 
for the finite lattice of A^ = 36 
sites shown in Figure [3l The 
site indices j = 1,3,4,6,9,13,24 
correspond to those of Figure [3l 



atoms — the simplest possible side groups. 

For haydeeite, the experimentally defined H position yields an unusually short O — H distance 
(d) of 0.78 A (for kapellasite the H position has not been reported). Therefore, in LDA 
calculations we relaxed the H position which resulted in an optimized d ~ 1 A for both systems. 
In order to elucidate the influence of d on the exchange integrals, we have studied the range 
d = 1.0 lb 0.3 A which covers the experimental errors (that are large due to weak x-ray scattering 
by H atoms). The details of our calculations are described in [3]. We find that d has a strong 
influence on the NN exchange Ji, while Jd is almost unaffected. Obviously, this modifies their 
ratio a and consequently the physical properties. The underlying reason is the shift of oxygen 
states down in energy caused by the shortening of an O — H bond. In consequence, Cu — O 
hybridization decreases. Therefore, for a fixed Cu — O geometry, a decreased d reduces the AFM 
contribution to the total exchange. This is supported by a strongly reduced bandwidth (which 
is related via tf^jf- to the AFM exchange) for a reduced d (Figure EJ right panel) and fits 



well the results of supercell calculations. There, we find a quasi-linear Ji{d) dependence with 
a positive slope for both compounds. Ji is 2.5 times larger for d = 1 A than for d = 0.8 A 
in kapellasite (2.5 meV versus 1 meV), while it changes the sign in haydeeite (0.7 meV versus 
-0.2 meV). 

The next step towards a physically relevant picture of real materials is to evaluate the influence 
of a onto the GS and the low-lying excitations. Here, we use a Ji-Jd spin-1/2 Heisenberg model 
on the kagome lattice {Ji-Jd model). The classical GS of the pure kagome system {Jd = 0) 
is known to be highly degenerate. The additional diagonal bond Jd reduces this degeneracy 
drastically and selects non-coplanar classical GS's with twelve magnetic sublattices (see Figure E]) 
among the huge number of classical GS's existing for J^ = 0. For the quantum spin-1/2 model 
we have calculated the GS by exact diagonalization for the finite lattice shown in Figure O 

Obviously, the quantum GS spin-spin correlations are drastically affected when including 
Jd- While for J^ = all spin correlations except the NN correlation function (SqSi) are close 
to zero we find a well pronounced short-range order for a > 0.25, that corresponds to the 
classical magnetic structure. For a > 0.55 the strongest spin correlations are (S0S3), (SoSe), 
(S0S9), and they belong to the chains formed by Jd bonds. The NN correlation function (SqSi) 
monotonously decreases with Jd- The correlation function (S0S24) belongs to two sites located 
on sublattices being perpendicular in the classical GS, and it is almost zero for a > 0.55. This 
leads to the conclusion that even in the quantum model the GS might have a non-coplanar 
magnetic structure giving rise to enhanced chiral correlations. 

To summarize, we have shown the relevance of the ligand field for the low-energy physics 
by example of two spin-1/2 cuprates with a kagome lattice geometry. We found that there are 
two relevant exchange integrals in both materials. The ratio a = Jd/ Ji is strongly dependent 
on the O — H distance and has a drastic impact on the physical properties: the quantum GS 
of a corresponding spin-1/2 Heisenberg Ji-Jd antiferromagnet on the kagome lattice exhibits 
strong magnetic correlations along Jd bonds for a > 0.5, whereas the other correlations remain 
weak. Therefore, the low-energy excitations might he S = 1/2 spinous causing an effectively 
one-dimensional low-temperature physics as has been discussed previously for other 2D models 
such as the crossed-chain model [10] as well as for the anisotropic triangular lattice [llj . As 
the physical properties are related to a which is strongly dependent on the O — H distance, the 
H position should be reinvestigated experimentally with better accuracy. Generally, a precise 
determination of the position of side group atoms is crucial for an accurate model derivation. 
We suggest that the dependence of exchange integrals on side group positions can be used to 
tune a magnetic ground state under changed conditions like external pressure. 
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